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CHAPTER - I 


IHTROPUCTIOI 


Water is a renewable source that follows in 
nature a path called the hyd rolo gic cycle. The details 
of the hydrologic cycle are described and analyzed at 
length in texts on hydrology( 1 , 2 , 3). We will limit our- 
selves here to point out that oceans are the sources 
of all waters on the earth. Energy originating in the 
sun and reaching this planet causes the evaporation of 
water from the oceans. Meteorological processes results 
in the precipitation of part of the atmospheric water 
vapours. Water then reaches the ground as liquid 
(rain, dew) or is solid state (snow) and eventually 
finds its way back to the ocean, thus closing the 
hydrologic cycle. 

The determination of the casual relationship 
between the application of water and its use, was the 
starting point of man’s attempt to top the hydrologic 
cycle at one or more points to utilize water for a 
variety of pruposes in his social practice. Man's 
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attempt to take advantage of certain aspects of the 
hydrologic cycle, gave rise to water resources projects 
of a wide range of sophistication. Ir their simplest 
form, such projects are nothing more than primitive 
facilities for the storage of rainfall in cisterns, for 
human consumption and domestic use. In the deltaic 
region where agriculture first developed, the river 
flood level being higher than the surrounding country, 
natural flooding for irrigation could be easily done. 

Open cuts were made in the river banks to bring, the water 
from lower levels and take it further away from the 
river. Once the necessity of control was realised, 
some primitive type of head regulator was evolved. 

When the rivers fell very low these open cuts could 
not draw any water, and then a sort of bund across the 
river made of wood and earth on stores was evolved from 
which sprang up the masonary weirs and dams. 


Modest advances in utilization of waters can 
be detected in the fourteenth century during the 
Sennaisance period. For the scientific study of water 
resources, we owe to Daniel Bernoulli, who in 1738 applied 
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the concept of oonseryation of matter to hydrodynamios . 

At about the seme time, d' alembert showed that Hewton's 
third law Of motion applies eiually well to moylng 
bodies (hence, to flowing water) as to bodies at rest. 
Toward the end of eighteenth century, Ohezy derived 
the formula describing the flow of water In open channels. 
An important laid mark in the development of modern 
thought in relation to water problems is the treatise 
published by Darcy in 1856 regarding flow of water in 
porous medium. A concise summaiy of the historical 
aevelopment of scientific thinking related to water 
resources can be found in Chow's, Hand-book of 
Applied Hydrology ( 4 ) . 


A common characteristic of most regions is a 
apace and time in balance of water demand and natural 
supplies. Han in his social practice finds that 
periods Of highest rainfall do not odBclde with periods 
of maicimum water demand, fable I illustrates the 
relative disparity between monthly water supply and 
irrigation demand in the Bari Doab Tract, Bunjab. 
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TABES I 

Montlily fluctuations in rainfall and y^ater 
demand in EBD tract, Punjab as percentage 
of yearly total (5). 


Month 

Rainfall 

Water demand 
for agricul- 
ture 

J an . 

5.6 

3.2 

Peb. 

6.1 

3.9 

March 

5.0 

7.5 

April 

2.8 

10.0 

May 

1.6 

13.1 

June 

6.9 

13.4 

July 

28.0 

11.5 

August 

29.6 

10. '8 

Sept. 

9.6 

11.0 

Oct. 

1.3 

7.5 

Nov. 

0.4 

4.1 

Dec. 

3 . 1 

3.8 

Total 

100.0 

100.0 
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Ho remove this time imbalance, man in his 
social practice has resorted to impounding excess 
water supplied during times of high runoff and 
releasing the stored water during times of demand. 

Historically, water projects were first conceived of 

as single purpose project. We have in the course of 
recent history seen the evolution from this concept 
to multi-purpose construction. The use of multiple 
purpose reservoirs in regional development schemes 
is today a well established practice (6)» 

During the last ten years, the development 
of water resources in any basin is being thought of 
as an integrated development (conjunctive use) of 
ground and surface waters. Consequently conjunctive 
use of ground water should be considered wherever 
major aquifers of good water yielding properties 
exist in a region in which large volumes of surplus 
water can be stored and withdrawn cyclicly. A 
distinctive advantage of storing water underground 
over a surface reservoir is the water stored is 
not subject to evaporation losses, furthermore 
savings may be achieved by leaving the canals 
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unlined, because seepage from canals percolates 
to the aq.uifer and is not lost. Such a dual system 
while certainly advantageous from the view of water 
conservation may require increased capital outlay 
for wells, reservoirs, pipe lines, canals, pumping 
station and artificial drainage facilities. A 
comprehensive survey of the pros and cons of 
conjuctive water use is given by Todd ( 3 ). 

A system under conjunctive water use may 
offer the alternatives of extracting water from a 
number of aquifers, by pumping, or from reservoirs 
by gravity flow and/or replenishing the aquifer with 
recharge facilities. Given all these sources and 
facilities the major question we will be faced with 
is to what extent each of them must be used on to 
what scale expanded to give the most desirable 
results. Modern operational techniques like linear 
programming and dynamic programming offer possibilities 
for finding solutions to problems of optimum timing, 
scale of development and use of alternative water 


sources. 
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The present study will use the problems of 

agricultural development ip Pup lab «=? t + 

j-u rup^ao as Its source of 

iPsplratiop lor the deyelopmeut of a determlplstlo 
linear Programming Model lor a water resouroes system 
eomprisiog of canals, tube wells and drainage and where 
irrigation is the only water aotlvlty. Water require- 
-nts Of the orops that are grown is met by a conjun- 
ctrye use of ground and surface waters thrcugh tube-wells 
and canals re spec tlyely. The model as developed will' 
include as input data, the hydrologic and economic 

characteristics of the water basin. Our objective 
Will be to determine the extent of allocation of 
inrlgated area to alternative orops and the seasonal 
water releases from the two sources, canals and tube- 
wslls, to meet the seasonal water requirements of the 
crops, for a one year period of operation. The policy 
we will be determining will be such that the returns 
function for the irrigation activity, subject to a set 
of constraints is maximised. Chapter III of this 
study deals with all these. 

In Chapter IT we will apply the model developed, 
to the Ban Doab Tract Comprising of Curdaspur and 
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Amritsar districts of the Punjab State. The Bari 
Boat Tract is situated between the two rivers Bavi 
and Beas of tlie Indus River system. ' The available 
hydrologic and economic deta for this region is 
sufficient for the limited scope within which a concep- 
tual model is developed in Chapter III. 

The results and conclusion of this study are 
presented in Chapter V. Appendix A gives the References 
and Appendix B the computer program used to solve the 
optimization problem posed. 



CHAPIEE - II 


EB7IEW OP EI1EBATUEB 

One of the first to state specifically the 
needs for conjunctive ground and surface waters use 
and suggest steps to he followed was H. Conklin ( 7 ) 
in 1948. He proposed the following conditions 

necessary to the implementation of conjunctive water 
use. 

G-round water of the alluvium is readily 
accessihle to extraction by pumps j 

The alluvium is sufficiently permeable to 
yield water in commercial quantities j 

The alluvium is recharged naturally with 
water or is susceptible of being recharge artifij*- 
ially to such an extent that heavy commercial 
drafts can be sustained? and 

If artificial recharging is necessary it can 
be done at feasible cost. 


4 ) 
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Conklin's work to help establish the fact not 
generally appreciated then, that "ground water Is 
merely water beneath the earth's surface and that 
eurfaoe water is merely water on the earth's surface. 
The same drop of water mowes from abowe the surface 
to below the surface and vice-versa. IShether above or 
below, it follows the same physical laws at all time", 
pen-d the doors for an integrated development of 
ground and surface waters. 

Ii6 efficient managecE nt of water resources 
continued to attract great attention in the face 
of increasing demands of water for variety of purposes 
coupled with the limited availability of natural water 
resource. An increasing necessity to consider 
water problems not at a local level but at a regional 
level and also an integrated approach for the utiliza- 
tion of the available water resources was felt. The 
problems as formulated continued to become more and 
more complex and available conventional methods proved 
incapable for their solution. Dantzig’s 
development of Linear Programming technique in I947 a 
and Bellman's (9) development of technique o.f 
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programr.iBg pro.lasa the mathematical tools 

necessary for the solution t. a.- 

^ vexing water problems. 

- use 01 system techniques and the advent ol 

computer ‘technology urovidp^ +v. 

, ,, Of impetus for a oonsl- 

erable extent/wori: towatds efficient wate 

iicient water management. 

Hall and Buras ( , 0 ) have outlined the 

potentials of om-i 

programming to multistage decision 
blems regarding allocation cf water to variohs 

alternative uses chnipp 

» hoice among alternative reservoir 

sites and optimum scale of 

^ 01 development of a water 

conservation system. 

'*at started as optimisation ol a lew partl- 
oular physical parameters, soon developed where 
economics played a major factor in the developm«.t 
Of water projects Involving optimum water resource 
allooatlons. Maass et al gr., , eomprehenslve 
fuformatlon of budgeting and investment policy and 
the relationship between income distribution and the 
economics of water resource projects. 

Buras in 1963 ( 11) used dynamic programming in 
peratio.i ol a three state system Including 
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one surface reserve i 


irs, one aquifer 


facilify and two 


9 one recharge 


irrrgatioh water would be 
reservoir or the aquifer. 


agricultural fields toward which 


diverted from tte surfj 


Fowler in 1964 (, 2 ) approached oonjunotiTe use 
problem IrOM the engineering point ol view. He sugeest- 
od the following steps toward the developa.ent of a 
operational plan for optimum conjunctive use of ground 

d SUTf3.C6 

1) A thorough geologic investigation has to be 

made and with the estimates of the geometry of 
the geologic formation and hydraulic aquifer 
properties an malog of the hydrologic basin 
can be constructed. 

2) The hydrologic history of the river basin has 

to be studied to understand the factors ielating 
to the hydrologic inputs and outputs. 

3) With the construction of an operational 
analog model of the system, an optimal oper- 
ational policy can be developed. 

optimum policy developed 
and to overcome legal and financial bottle 
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neckB, a,: oparatiag agapoy with adeguata powera to 
oonttol or cooperate in the control of the warloua 
available water reaouroes has to be established. 

Chun at al in 1964 03 } aescribe the develcp- 
-Pt ot an analog aodel of the water supply, .latn- 
hution and replenishment system of the los Ungeles 
^-in. Ihe basis was divided in to 81 polygons of 

varying sizes and the parameters of the hydrologic 

equation were fpfi in +r. ^ 

fed inoo a computer for each of the 

polygons. Using historical data collected over several 

yoeuo, the simulated behaviour of the analog was tested 

against the historio behaviour of the physical system. 

Ihls work Is probably the most comprehensive and 

detailed regional water plan. 

Wille considerable work was being done in the 
oonjunotlve utilisation of ground and surface waters, 

legal aspects assumed importance. Valentine in 1964(14) 

elves a lucid commentary on the legal aspects, const- 
raints ana approaches to be considered in the planning 
Of conjunctive water use. According to Valentine 
the uncertainties in the knowledge of hydrologic Processes 
10 the major reason for the lack of a comprehensive 
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ground water law behind the establlahed surface 
water legislation. 

The los Angeles basin study no e 

v15; includes an ou1 

standing economic analvsiq Wom 

halysis. However due to the relative 

unimportance of aa-rioni+,, 

agrieulture xn the Lavy Vegas Valley, 

=tuay did not include the consideration of the 

-.Pier problems of egrloultural flexibilities. 

Considerable worh has been done in the area of 
conjunotive use and for situations analogous to the 
scope of the present study, by the Center of Popul- 
ation Studies, Harvard university. Pouglas V S.ith 

1970 (16) giTes a lucid and comprehensive report 

on irrigation Planning Models anr^ -t-h 

6 mu a 0 IS and the present study 

owes considerably to this woifc. 








CHAPTER III 


COHCBPTUAl' MOREL 

5 • T ' ihtrol tjo tioit 

In this chapter we will develop a conceptual 

model for the conjunctive use of ground and surface 

waters for a single purpose viz., irrigation. We 

will consider a one year period of operation consisting 

of monthly seasonal intervals. Our objective will be 

to determine the extent of acreages for the various 

crops and the seasonal water releases from the two 

sources, canals and tube wells to meet the seasonal 
water req.uiremen ts of the crops. The policy we will 

he determining will he such that the returns from 

irrigation subject to a set of constraints will be 

maximum. In developing the model we will make an 

assumption that all input data are deterministic in 

nature. Pigure 3.1 gives the schematic representation 

of the model we will be developing. 

3.2 SYS TIM 0 OMPQCT TS 
3.2.1 Cropped Area; 


In the cropped area different crops having 
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(Jifferent seasonal water requirements, different 
yield and different value per pound are grown. The 
cropped area is subdivided into two zones, irrigated 
and non-irrigated areas. Natural precipitation in 
the area is available to meet part of the seasonal 
water requirements of the crops. Additional water over 
this amount is to be met by the seasonal water releases 
from canals and tube wells. Of the water delivered to 
this area, part of it will be lost as surface run-off, 
aquifer recharge and non— bene fi cial svapo transpiration . 

3.2.2 0 an al s ; 

Ihe canal system takes off from the river at the 
canal headworks and delivers water to the irrigated 
area. In the period of transit from canal head works 
to the irrigated area, part of the water that is diverted 
to the canals will be lost as surface runoff, aquifer 
recharge and non-beneficial e vapo transpiration . 
Engineering and economic factors vifill impose a 
maximum value lor the canal capacity. The canal 
system is not bound to be operated during the entire 
season. Due to routine and also unexpected necessities 
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of repairs etc», the efficiency of the canal system 
will be less than unity, luring a season the demand 
from canals will have a peak period and hence the 
ratio of peak to average demand is also to be consi- 
dered. The operations cost lor canals is given in 
rupees per acre feet of water diverted to the canal 
system. 

3.2.3 Tube-we lls; 

The tube well system pumps water from the 
aquifer and delivers it to the irrigated area. 
will also consider the possibility of the tube well 
system pumping water direct to the drainage system. 
Similar as in the case of the canals, the tube well 
system too will have a maximum capacity imposed by 
engineering and economic factors, will not operate 
during the entire season, will have an efficiency 
less than unity and the ratio of peak to average demand 
will also have to considered. Although strictly the \ 
operations cost for pumpage will be related to both | 

(Juantity and extent of lift, in our model we will | 

assume this cost as given only in rupees per acre-feet 
of water pumped. 



The drainage system collects the surface 
runoff from the canal zone, irrigated and non-irrigated 
areas and the direct pumpages into it if any from 
the tube wells system and delivers it to the river. 
Similar as in the case of canals and tube wells, the 
drainage system too will have a maximum capacity 
imposed by engineering and economic factors, will 
not operate during the entire season, will have an 
efficiency less than unity and the ratio of peak to 
average demand will also have to be considered. 

Yftiile delivering the water to the river; as pumpage will 
be necessary, the operations cost for the drainage 
system will be strictly related to both quantity and 
extent of lift. However in our model we will assume, 
that this cost is given only in rupees per acre-feet 
of water delivered. 

3o2‘5 A quifer s 

It is feasible to exploit the underlying 
aquifer system for irrigation use by pumpage of water 
from it through a tube well system. Geologic, hydrologic 
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c.iiu ecciioiiic factors will impose a maximum value for 
one annual pumpage from the aquifer system. Natural 
aquirer recharge from adjoining aquifers, aquifer 
j-uchouge from canal zone, irrigated and non-irrigat ed 
areas form the inixit into the aquifer system. Besides 
the water that is pumped out by the tube well system, 
\v..o.;.r !,'ill be also leaving the aquifer in the form of 
evaporation and aquifer return flow to the river. 


SYST EM PARiil EllBB : 


Having explained the system components in the 


preceding section, V7e will now proceed in this section 
to define the system parameters. 

3 «3 .1 P.=‘^,9PPed..A^^: 

-'’•-''*1 ‘I . Input Variables: 


(i) TC 

(ii) TIA 

(iii) Xj 


- total cropped area in acres 

= total area available for irriga.tion 
in acres 

= fraction of . bhe irrigated area 

allotted for the jth crop. j= 1 , . . .J 
v?here J is the total number of crops 


grov'/n. 
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(iv) 



(v) YPA. ■ 

J 

(vi) VPP 

«J 


water requirement of the jth crop 
during the ith season, ih inches, 
i= 1 , ♦ . . I, where I is the total 
number of seasons and j = 1 , . , , J. 

yield per acre of thf- jth crop, in 
pounds per acre, j = 1 , * ■ . J 

value per pound of fche jth crop, in 
rupees per pound, j = 1 , , . . J 


(vii) PI^ = natural precipitation in the area 

during the ith season, in inches. 

(viii) SRg, AR2 > ET2 = fractions of y/ater delivered to the 

irrigated area that is lost as 
surface runoff, aquifer recharge 
and non beneficial evapotranspiration 
respectively - 


(ix) SR^, idilj , Elj = fraction of water delivered 

to the non-irri gated area that 

is lost as surface runoff, aquifer 
recharge and non-benef icial evapo- 
transpiration respectively. 

(x) ©2 “ fra.ction of water delivered to 

the irrigated area that is used to 
meet the water requirements of 
crops .02=1 -SR2“AR2-ST2 • 
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3*3.1 .2 Operating Tariables to be Determined; 


(i) 

^3 = 

irrigated area of the jth crop 

lr3..2 

Canals; 


3.3*2 

.1 Input Variables; 

(i) 

Xi 

discharge in the river at the canal head' 
works during the ith season in acre-feet, 
i = 1, I. 

(ii) 

Y = 

maximum capacity of canal system, in 

cusecs. 

(iii ) 

II 

s 

number of days the canal is operated in 
a season. 

(iv) 

II 

efficiency of the canal system 

(v) 

^(y) , 

ratio of peak to average demand of 
the canal system. 

(vi) 

o 

II 

operant ions cost for the canal system 

in rupees per acre-feet. 

(vii ) 


= fractions of water diverted to the 


canals that is lost as surface 
runoff, aquifer recharge and non~ 
beneficial evapotranspiration respective! 
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(viii) e, = fraction of water diverted to canals 

that is delivered to the irrigated area 
= 1 “SR^ . 

3»3*2^2 Operating Variables to be Determined" 


(i) 

^i 

- flow diverted to canals during the tth 



season, in acre-feet, i=1 , . . . I. . 

3.3»3 

Tube-wells: 

3.3.3 

-1 Input 

Variables: 

(i) 

T 

= maximum capacity of tub ew ell system, 



in cusecs. 

(ii) 

a(t) 

= number of days the tube-wells are 



operated in a season. 

(iii ) 

^(t) 

= efficiency of the tube well system 

(iv) 

^(t) . 

= I'atio of peak to average demand for 



tube wells 

(v) 


= operations cost for tube v/ells in 



rupees per acre-feet. 
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5 *2 *3 *2 OpBraiiing '/axiablsg "fco b© DeljBmixisd. s 


(i) 

ti 

pumpage tnrough tube well system for 

irrigation during the ith season, in 

acre-feet, i = i ^ 

(ii) 


pumpage through tube well system direct 

into the drainage system during the ith 

season, in acre-feet, i = 1, I, 

3.3.4 

Drainage ; 


3 . p • 4 • 

1 Input Variables; 

(i) 

Z 

maximum capacity of drainage system. 



in cusecs. 

(ii) 


number of days the drainage system is 



operated in a season. 

(iii) 

= 

efficiency of the drainage system 

(iv) 

= 

ratio of peak to average demand of 

drainage system. 

(v) 

= 

operations cost for the drainage system 


IS rupees per acre-feet. 
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5»3» 5 Aquifer; 

3«3*5.1 Input Variables; 


(i) AMI. 

(ii) AEP 

(iii) iili’.E 

( iv ) pm 


= annual aquifer natiiral recharge from 
adjoining aquifers, in acre-feet, 

= annual aquifer eTaporation, in acre-feet 

= annual aquifer return flo\? to the rivex’, 
in acre— feet 

= permissibUe mining rate in acre-feet per 


year , 


3,4 miHSMATIG AL PO E?U71I IOF: 

Having explained the system components and 
defined the system parameters in the preceeding sections, 
v/e ¥;ill in this section proceed to develop the mathematic 
f ormli tion in the form of a Linear Programming Model. 

3,4,'l Ob j ecti ve F unction; 

The objective function includes returns from 
the irrigated area and the operations' costs for canals, 
tube wells and drainage 

3.4.1 ,1 . Return Prom Irrigation; 

Annual returns from the irrigated area is the 
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value of the 


crops giovn there. This will be equal 


to 


J 

z : 

0=1 


Au z ypi X 7PP. 

J J -J 


»2 Operations. Cost P 


Iror Canals: 


Annual operations costs will he equal to 


I 


H 

i=1 


(y) 




3.4.1 .3 Operations Cost For lube-wells; 

Annual operations costs will be equal to 


(t ) 


i=1 


z (t^ + s ) 


3. 4. 1.4 Operations Cost Por Drainage; 

Annual operations costs will be equal to 


i=1 


(z) 

C X z . 

1 
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wnere z. is the water delivered by the drainage 
sj-steni during the ith season. 


+ SR. 

1 3 C Ji + 

sRg (e^ 

yi + 


+ (PIj_ X 

SE2/I2.O) 

J 





0=1 



+ (H^ X 

SS3/I2.O) 

(PC - 

J 

H 

^ 3 > 




3=1 


' ?! (.SB., 

0. SE 

) + t^ 

SRg 

+ s. 
1 


J 




+ (PIj/12 

. 0 ) SI 
3=1 

A. (SR, 
J c 

p - SR 

3) 


+ (Pl^/ 12 . 0 ) X TC 

V/e see that 2 . is represented as a function of 
^i? tj^s and A. .. 

The objective function to be maximised nav 

becomes 


J 

A . X YPA. X 

0=1 


- n {o(y> y, + 


i=1 


(ti+sy 


+ 0*-^^ X z^) 
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■l«A»2__Co_nst ralnt Equa-oionss 

Tae above objective function is to be maximised 
subject to the following constraints. 

fii® seasons.1 diversions to canals cannot 
exceed the seasonal inflow into the river at 
the canr.l hesd v^orks* 




ii) 


v.ater divei'ted to canals in any season cannot 
exceed the canal capacity 

S) X X j^/ X Y 

where -9 is the conversion factor from acre-feet 
per day to cusecs. 

iii; iiciber punipeo. from tube wells in any season 
cannob exceed the tube well capacity 


i?x 


(t ) 


( Si + "fei ) / d 


(t) 


X 




(t) / 


viv) ater diverted to draiiiage system in any season 

cannot exceed the dPctincige capacity 

X X z,. / d*^^) X Z 
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(v) Total area of Tarious crops cannot exceed 
“the 1^01:81 available area for irrijgaiion* 

^ A. < TIA 
d=i 

(vi; Total water removed from aquifer in a year 
cannot exceed the mining rate 


i=1 


(s^+t^) - AE^x y^ - AR^ X (6^ x y^ + 


J 


~ (PI.:/12.0) X AR^ ^ ^2 ~ '^-^i/''2*0) X im^x(TC-.^ 

d=1 

+ AST + AEP - AAE ^PI,IR 


(vii) Water requirement for the crops is met in 
each season. 


J 

2S h ^ "Wj - ®2 (©1 yj. + tp 

,1=1 


J 

“ 51 ^2 ^ ®2 ^ Plj_/12.0 

d=T 



0 


(viii) Hanagement considerations restricts a maximum. 

value for the irrigated acreages under each crop. 



• ita Cl, conceptual model for conjunctive 
use of ground and surface waters for a single purPOse 
vxz., irrigction thus developed, we will proceed to 
Ux-plj it to a physical model in the next chapter. 






Jammu and Kashmir 


jurdaspur/, 


Amritsar 


N Hoshiarpur 


JuUundur 




/ Feroze.pur 


FiG-4d .INDEX plan OF BARI 


chzipt.'eir - rv 


PPIYSICAL I'iODEL 


4.1 IFTROroC'JlOij; 


Ban Boab Bract in Punjab State of the Indian 
Union forms the physical mcdel of our study. This area 
includes the tivo districts of Amritsar and Gurdaspur. 

She area is bounded on almost three sides by the Hart and 
Beas riyers of the Indus rlyers basin and by lakistan 
on the fourth - western side, ligure 4.1 gires the index 
plan o£ the tract. Topographically Gurdaspur district 


lies in the nerrow strip of territory stretching between 
me Himalayas and the Indo-Gangetic plains into which 
run the spurs of the Himalayas and Aamitsar district lies 
in me Inao-Gangetlc plains. The height of the Doab in 
the e.rtreme Iforth-Hhst is 1200 feet abore sea lewel and 


ill the Southwest 


is 685 feet and slope from Borth-East 


to South-17est, 


Ihe totaa area of the tract is about 3,300 
sq .miles. The population as per 1971 census was 30.45 
lakns. Out of this the rural population was 22.65 lakhs 
and urban population 7.80 lakhs signifying the dominance 
of the agricultural sector in this area. 
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4*2 I IBTaOEOLOG-I CA L GQlT n7T^TnT7S : 

4*2.1 Climates 

Due to the competative paroximitj- of the hills 
and the prevalence of canal irrigation, the area is 
temperate in the hot months from May to September. There- 
alter the air becomes drier and cooler every day. Prost 
is comm.on during January and February. The temperature 
r. nges from a la? of 40*^F in the winter months to a 
high of 1 1 0°P in the summer months. 

i^2_.2_Rai_itf a^lj_ 

The general dis cribution of rainfall over the 
area is given infigure 4.2. We notice that rainfall 
decreases from Srorth East to South West. Rainfall 
figures for Gurdaspur and Amritsar districts for the 
years 1902-03 and 1962-65 are given in Figure 4.3* 

The average annual rainfall for Gurdaspur district 
works out to 37*10 inches and that for Amritsar works 
out to 22,72 inches. Figure 4*4 gives the distribution 
and percentage of total annual rainfall over the 
months for Amritsar District. 
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4.3 S EOHYDSOLGCtY; 

The hoab 1ms been formed, by alluvial deposits 
brought by the rivers iiavi and Beas. Separate GeO“ 
hj'drologic studies for G-urdaspur and Amritsar Districts 
have been carried out b;' the Geological survey of India 
c-nd extensively reported by Punjab Government 
publications (18) and by Uppal et , , al(l 7 ). These 
studies reveal that the substrata mostly contains medium 
sand . There is generally 15 to 20 feet deep soil 
crust at the top and is follovYed by water bearing strata 
CO. misting of medium or coarse sand mixed with very 
fine traces of clay. 

Groiindwater in the area occurs under confined as 
we].l cm iinconfined conditions. ’later table is at the 
lowest level in June before the mmnsoon season and is 
at the highest level in October after the monsoon 
season. The principal sources of gTOu d water acre 
infiltration from rainfall, seepage from irrigation 
channels and irrigated fields and the subsoil flow from 
upper parts. Ground vjater in the area is generally 
fresh and Bti.idies (18) indicate that large scale 
d.evelopmerit of ground watei' through tubewells is feasible 
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in a grei.t part of the area. Table 3.1 gives the 
storage capacity of the area. Table 3.1 gives the 
store ge capacity of the ground water reservoir as 
reported by Uppal et . , al. 


TABES 3.1 

STOIlAGrS CAPACITY OP G-ROITID '-'fATER RESERVOIR 



IE THE 

DOAB TRACT 


jjepth gone 
(in feet) 

Thickness 
of strata 
(in feet ) 

Average 

Specific 

Yield 

(fO 

Volume of 
wat er 

(Ilillion acre 
feet) 

15 ^ 50 

35 

12.9 ‘ 

9.5 

50 - 100 

.50 

17.5 

18.4 

100 - 350 

250 

17,5 

91 .9 
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In the conceptual model we had developed in the 
III Chapter, the aquifer system requires as input data 
the followings annual evaporation, annual sub-surface 
return flow to the river, annual aquifer recharge from 
adjoining rivers and a maximum permissible mining rate. 

All these parameters appear in the mathematical model, 
in the constraint equation specifying restriction on 
the yearly mining rate (refer 3.4.2). The same 
constraint when expressed in the standard manner has the 
right hand side constemt ass 

PMR - ARP -..ART + AAl 

J 

+ (PI./12.0) X AS^ X (TO - A.) 

^ . pi ^ 

■\?hile applying the mathematical model to the physical 
model we will assume that, ’’PBffl - ARP - AST + AAR” = 5.4 
million acre-feet. This is a reasonable estimate consi- 
dering the ground water capacities for different depths 
given is Table 3.1. A comprehensive water balance of 
the area would lead to a refinement for this value in 
future studies. 

AGRIOUITURE s 


Iheat, maize, rice, sugar-cane and cotton are the 
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main crops grown in this area. The present level of 
cropped area for each of the crops is given in Table 4.2. 
Besides these main crops^ we have included all other 
crops in one single head. Table 4.2 also gives the 
percentage of the total cropped area for each of these 
crops. While developing the constraints on the maximum 
permissible irrigated area for each of the crops, we will 
mahe an assumption that the percentage of the total area 
available for irrigation for each crop would be the same 
percentage value given in Table 4.2. 

TABEB 4.2 

CHOP CROPPSD AREA PRACTIOH OP TOTAL 

OP EACH CROP CROPPEP AREA 

(in million 

acres)., 


Wheat 

0.96 

0.40 

Rice 

0.38 

0. 16 

Maize 

0.24 

0. 10 

Sugar-cane 

o 

0 

V 

0.06 

Co tton 

0.08 

0.03 

Other Crops 

0.60 

0.25 


Total 


2.40 


1.0 
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CHAPTER Y 


RESULTS AND CONCLUSION 

5.1 GEHSEAL? 

Tile objective function of the Linear Programming 
Model we formulated had 42 operating variables to be 
determined; the extent of irrigated areas for the 6 
Crops we' considered, the monthly releases from tubewells and 
canals for irrigation, in a year of operation and the 
monthly releases from tube-wells direct into the drainage 
system. The number of constraint equations for the 
Linear Programming Model we formulated was 68; 12 
constraints composed by restricted river discharges, 12 
constraints imposed by restricted canal capacity, 12 
constraints imposed by restricted tube-well capacity, 12 
constraints imposed by restricted drainage capacity, 12 
constraints imposed by the necessity to meet the water 
requirements of the crops, in all the .above cases one 
each for the 12 months in a year's operation, 1 const- 
raint imposed by restricted mining rate from the aquifer, 

6 constr ainig impo sed by restricted area available for 
irrigation for the individual crops considered, and 1 
constraint imposed by restricted area available for 
irrigation. 
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Simplex alogrithm was used to determine the 
optimal policy. Each run of the Linear Program where 
the Simplex tableau was generated as a single 
array matrix from the required hydrologic, geohydrologic 
and economic data fed, took about 4 minutes of 
IBM/7044/1401 computer time to reach the optimal 
solution. 

5.2 OPTIMAL POLICIES ; 

Ihe results of the optimal policy arrived at 
are given in Figures 5.1 and 5 . 2 . 

5.2.1 Opt i mal Water Releases s 

is would be expected from the prevailing 
comparative costs of canal and tube-well waters, 
we find from figure 45.1 that canals are used to 
maximum capacity of 202 thousand acre feet in all the 
months where water release for irrigation is required. 

We also find that tube-wells are used to full capacity 
of 404 thousand acre-feet in the month of October and 
in all other months when operated, at varying amounts, 
all less than this value. In the monsoon months of 
July, August and September we find that all the water 
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requirements are met by natural precipitation alone. 
Figure 5.1 also gives the net water made available 
for irrigation after abstraction of all forms of 
losses considered in the model. 

5 . 2.2 Optimal Cropping Policy; 

Figure 5.2 gives the optimal acreages of irrigated 
area under each of the 6 crops we considered in the 
model. Y/hile wheat, maize, sugar-cane and other crops 
enter the solution at their maximum specified acreages, 
rice centers the solution at a level of 0.215 million 
acres, which is less than the maximum specified level 
of 0.288 million acres and cotton does not enter the 
solution at all. 

5.2.2 Total Benefits^ 

The total benefits accrued from the irrigation 
activity as per the optimal plan is 982 million 
rupees. 

5 .5 S31TSITIVITY MALYSIS ; 

The linear programming model also provides 
information about the sensitivity of the results to 






































errors in estimation of toe input data. Sensitivity 
analysis were carried out on the tube-well capacity, 
on the area available for irrigation, on the operation 
costs for canals and tube wells and on the values per 
pound of the crops considered. .Results of the 

sensitivity analysis are reported in Figures 5.5 to 
5.6. 


Cana l Oapacitys 

^ig^re 5.3 reports the changes in the optimal 
ution over a range of canal capacities. For all 
the^ cases considered we find that there is no change 
in t e canal releases. A significant increase in . 
the tubewell release with increase in tube well 
capacity is also found only in the months of October 
and November, meat, maize, sugar-cane and other 
crops enter the solution at their maximum specified 
level in the model at a caaal capacity of 12,000 
cusecs and with increased capacity of tubewells 
thereafter, they show no change. Increased canal 
capacity from 12,000 cusecs onwg^ds only forces 

irrigated area under rice to enter the 
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solution. At a tubewell capacity of .18,000 cusecs ' 

rice enters the solution at its maximum specified 

level of 0.288 million acres, and cotton also enters 

the solution at a level of 0.014 million acres. I 

At a canal capacity of 20,000 cusecs, cotton also 

enters the solution at its full specified level and 

the entire area of 1.8 million acres available lor 

irrigation is irrigated. As can be expected, the 

benefits increase with increase in tubewell capacity 

and are all shown in Pigure 5.5 The marginal benefit 

for an increased capacity of 2000 cusecs is steady at 

28 million rupees in the range considered. 

5.3.2 Sensitivity A nalysi s on Ar ea A vail^le Por 

I rri g ati on i;_ • 

Pig. 5.4 reports the. changes in optimal 
solution over a range of area available for irrigat- 
ion . Por all the cases considered we find that 
there is no change in the canal releases. With 
increased areas available for irrigation, the tube- 
well releases exhibit an increase in level and at- 
a level of 2.4 million acres, incidently the total 
cropped area in the region, the tube-well is used to 
full capacity also in the month of November. In all 

C£MT!IAL UBriA^IY I 

iu ^ 
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the oasee considered we find that maize, sugar-cane 
and other crops enter the solution at their maximum 
lamel speoilied in the model. In the case of wheat 
it enters the solution at its maximum specified level 
upto 2.2 million acres being made available for 
irrigation. At a level of 2.4 million acres there is 
a drop in the irrigated area under wheat as compared 
to the level of 2.2 million acres. In the case of rice 
there la a progressive reduction in irrigated area 
with increase in area available for irrigation. In all 
the oases considered cotton does not enter the solution 
at all. There Is an Increase in the benefits from 
irrigation activity with increase in area available for 
irrigation Implying an under utilization of available 
waters at lower levels of areas made available for 
irrigation. 

on Or)erptlnr>od Costs T^nr 
Q an al a n d 0?ubewells 


Figure 5.5 reports the changes in optimal 
solution for a range of operations costs for canals 
and tubewells keeping all other parameters same. 
For operations costs for canals at Js. 10 and &. I5 
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per a.cre-feet and for all “tlie operations costs for 
tube— wells considered we find that there is no change 
in the water release policy. When the canal operations 
costs is fis. 20 per acre-feet we find that the releases 
from canal enter the solution only in the months of 
October in its full capacity and in November, lor the 
same case, we find a marked increase in the tube-well 
releases during all the months, lor all the cases 
considered there is no change in the optimal cropping 
policy. As can be expected, the benefits show a 
decline with increased operations costs. 

5 >3 >4 Sensitivity Analysis lor Value of Props; 

ligure 5.6 reports the changes in optimal * 
solution for a range of values of individual crops with 
all o'ther parameters remaining the same, lor all the 
cases considered for wheat, maize rice and sugar cane 
and other crops we find that there is no change in ihe. 
water releases policy, lor the case when the 
value of cotton is Hs. 1.00 per pound we find that 
there is a difference but not a significant one 
for the tube— well releases while there is no change 
in the canal release, lor the same case we also find 




y:zi ;«y.i. .:ai 


AU^ALL 


AlAl C'A: ^CriT ,a.:l)rCKi) TC 


'III'- FOR^l'il 


SSASOBAl m* LEASES FRi^: CAHALS ' A® TUmSLI^ IR iTiOl©AM AI' iBa 

OF^rOI. ^LAK \ 


Other Cjx^g 


I-OLICY T!t FTLLICF ACR' 


iT'n' 


25J- - g : : :■ 

























44 


that cotton enters the solution at its maximum, 
specified level of 0.054 million acres and this is - 
at the cost of ric 3 whose irrigated area drops to 
0,205 million acres, For all other cases there is no 
change in the cropping policy. As cm be expected, 
with increase in value of crops, there is an increase 
in the total benefits. 


3.4 CASS Of RAIEFALL BBI1& QgB POUR’IH OP NOmiAL RAIKMI L; 

A case where the rainfall is assumed to be 
one fourth of the normal rainfall used in the model was 
also considered. Sor this case it w'as also assumed 
that the canal can be used only up to half its capacity 
we had considered in the model. The results of 
optimal solution for this case is reported in Sig. 5*7. 
¥e find that the canal is used to available capacity of 
10 thousand acre-feet in all the months of the year. 

Of significanee is the appearance of ceinal releases man- 
soon months of July, August and September. Tube-well 
releases also appear in all the months and the maximum 
capacity of 404 thousand acre-feet is used in the months 
of January, June and October. Sor all other months 
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tube well releases are lound to be greater than In 
the original optimal plan. Wiile irrigated areas of 
maize, sugar cane and other crops show no change in 
this case as compared to the original plan, we find 
that for this case, -there is a reduction in the areas 
under irrigation of wheet and rice. As can he expected, 
there is also a reduction in the total benefits for this 
case as compared to the original optimal plan. 

5,5 C PNG lUSI PITS ALTD SCOPE! 01 ? PUEIHBR STUhlSS 

Por a set of deterministic input data we had 
formulated a linear Programming Model which gives the 
water releases from two sources, canals and tube wells, 
to meet the water requirement of crops and the extent 
of irrigated acreages under each of the crop considered, 
so that the returns from irrigation activity is 
maximized subject to a set of constraints. This model 
was applied to the Sarl.Doab Tract in Punjab and the 
optimal solutions have been reported in the proceeding 
sections of this chapter. 

’That strikes one at the first glance of the 
optimal solutions presented is the absence of releases 
from tube— wells direct in to drainage system although 


4-6 


this was included in the model. 


One can easily explain 


It by the presence in the model ot only a cost factor 
lor this release and no benefit factor. ®ille 
planning the optimal policy lor a broader basin, say 
the whole of Junjab State, a benefit factor will get 
introduced for these releases since they can be made 


use of for 


irrigation in the lower regions. 


This 


would force their entry into the optimal solution. 

Another method, not considered in the model we had 
developed, to force their entry into the optimal solution 
would be to impose a minimum constraint on the amount 
of water to be pumped out of the aquifer in a year of 
operation. This would call for a comprehensive data 


on the ground water balance. 


In the sensitivity analysis carried out for the 
area available for irrigation, we found that with 
increase in area available, there is a reduction of 
irrigated area for rice and wheat in the optimal plan. 
Again this was mors pronounced in the case of rice. 

This is because the model we formulated does not 
include a minimum constraint for the irrigated area 
of individual crops. . Yfith the current accent on 
increased food grains production, it would be necessary 
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to include in the model minimum constraints for the 
irrigated acreages of these two crops.. This would 
ho?/ever call for a comprehensive demand function for 
these crops. 

None can fail to notice that cotton does not 
enter the solution of the Optimal plan, the reason 
being absence of minimum constraint for irrigated areas 
of different crops in the model we formulated. This 
however does not assume much importance because if the 
past trend is any indication there has been a progre- 
ssive reduction in the cropped area of cotton in the 
region of our study. 

The model we formulated had all input data as 
deterministic in nature, this^not truely reflecting 
the stochastic character of river inflows and natural 
precipitation, lor more realistic picture it would 
become necessary in future studies to incorporate the 
stochastic character of these variables. 

Another extension would be to include a 
constraint on the quality impairment. This constraint 
would hold the annual quality detirioration of ground 
water to a certain level. 
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Another extension would be to include a 
surface reservoir upstream of the canal head work. 

For the region of our study this assumes importance with 
the Ihien Dam coming up across the Ravi river in 
Himachal Pradesh. 
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computer program 


SIBJOB 

SIBFTC 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

G 

C 

G 

C 

G 

G 

G 

G 

G 


ir^HE main PPOGrAm THE SIMPLEX TABLEAU IS GEnErATED AS 

i, 0 ’?Sk%p?iMlL” 3 ;LGn 0 N^Is" 0 BfAIuED I„ SAME STOP AGE L oCAT lON 

N IS number of VARIABLES To BE DETERMINED IN THE OBJECTIVE 

FUNCTIONy^BER OP CONSTRAINT EQUATIONS OF THE. TYPE LESS 

than or EQAL to 0 , 1 M II 

ZS = 1 SS=M1 + 1 IPRl = N IPR2=M1 JL^-Ml 

SRluj)=iATkk'"RESuTRkMiNT’lN ItlcSs OF Crop J IN SEASON I 
YLDpA ( J ) sY IELDI N POUNOS PBR AGpE OF GROP J 
VAIPPC J)=VAlUE in rupees pEr pound of Grop j 
CAPC^CA pT.CAPD ArE CAPACITIES OE CAnAl .TUBE-WELL .AND 

OCclSa^AND Ocl^ARE OPERATIONS COST OF CANALS, TUBE V ILLS, 

°A^“"Ef ARE^FrcnONr^F^W^ATirkosT AS SURFACE RUNOFF . AQUIFER 

aIt^^AAr’ ARF "ArE'ANNuirAku!FfR°E«prRAnor annual AQUIFER 
RlcHApfE’F^^M AD^JoInUG koS?FERS^ AND ANNUAL AOUIFER RETURN 

PLOW RESpEglVELY qF GAnAi S, TUBE-WEllS AnD DrAInAGE 

RpT;ocr;pTAST""^ RPTAdS ArE RATIo o'f pEA< To AVArAGE DFmAND oF 

k«^ct DA^s?^D^YSD’ARf NSMBiroF DAYS CArAlS. TUBE-WEllS, 

p"„R ?f pk^MlSS^lrMlNlNG R^Tk ^E-FEET PER YEAR FROM AOUIFER 
TOTA is total cropped area TY-AvTr,M 

ttA t TnTA! area AVAIl^SlE FqP IRRIGATION tz-a-tt 

^RArUI ^^F^ACTI 0N ok ToTAl ArEA AVAIlABuE For irrigation 

made available for crop j .r c ccpt pwFr, DAY Tn Gl SECS 

Anew is ConVEpSIoN factor Fqr ^CrE-^Et pEr day to CISLCS 

DIMENSIONPRI hOO) »PR2tlD0) J/ f pp , - , ,Yi DpA ( z ) ,UP I A ( i p ) 

DIMENSION rIVroF(12)»PI‘12)»Wr{6M2)»VALPP(a)»YLUP 
DIMENSION PErAR ( 6l 

COMMONA(3500) ca. . 

Integer zs,ss,v,r»s,pri»pr2»fall 

real max 



R E AD 5 0 ’ ^ S s S S » M » M 1 
50 FORviaT (415) 

READ 6Q » I PP1 > I PR2 » JLl > JL2 
60 FORMAT (415 ) 

PEAD70 »MnS,MOC 

70 FormaT(2I5) 

read 71 » (pIVroF ( IA) , IA=i , wnS) 
read 71>(pi(IA)»IA=1? NrO ) ■ 

READ? I » 5 ( rtp ( J ojKO) »KO=l »NOS ) » jO=l sN0C3 
READ 71 » ( YlDPA ( JP ) ) JP=1 » NOC ) 

READ71 ’ ( VALPP ( JP ) » JP=1 »NoC ) 

READ 71 >CApC >CA pT ^CApD,cCC »nCT »0CD 
READ?! » t pFraR ( I) ? I = 1 »NnC ) 

?! Format ( 6Fio. 2 ) 

READ 72>SRl ,AP1 ,ETl 
read 72,SR7»AR2»ET2 
read 72»Sp3 ,Ap3 »FT3 
PEAD 72»AET»AaR >ARF 
read 72>EFC,EfT,EFD 
READ 72»RPTADC»RPTADT»rPTADD 
read 72»DAYSC,DAYST»DAySD 

72 FORMAT{3F10,2) 

READ 73»PMP 

73 FORMAT ( FiO. 1 ) 
rEAD74>TIA,T0TA‘ 

74 F0RMAT(2F10.1) 

ANEW=4840.0*R.O/( 24.0*60. 0*60.0 > 
CONC=RPTADC*AmEW/ (DAYSC*EFC ) 
ConT=rpTaDT*AmEW/ (DAYSTi^FFT ) 

CoND^RPT ADDijAmEW/ ( DAYSD*EFD ) 
theta 1 = 1 .0_SRi -A ri-ETI 
THETA2=1 .0-SRr-AR7-ET2 
C0ND1=C0ND^{Sp1+SR2*TheTa1 ) 

COND2 = COND^fSR2 

TAPI=0.0 

D075I=l>N0S 

75 TAPI=TApi+Dl(i)/i2.0 
A( 1) =0.0 

D0131I=1 ,NOS 
IR1=I+1 

131 A( Ipi ) = Rl\/onF( I ) ■}!■ tO.C i:-* 6 
D0132 1=1 »N0S 

IR1=IR1+1 

132 A(IP1)=CAPC 
D0133 1=1 »N0S 

r R 1= I Rl + 1 . 

133 A ( IRl )=CAPT 



D0134I=1»N0S 

IR1=IR1+1 

134 ACIrD = card - COND*Sr3*ToTA»PI ( n /i2.0 
IRlsIRl+l 

A(Ir1)=Tia 

IR1=IR1+1 

A< I R1 ) =Pf^R+AAp-AET-ARF+ToTA*TAP I*AR3 

D0135I=1 »NOS 

IR1=IR1+1 

135 A(IR1)=0.G 
D0136I=l »NOC 
IR1=IR1 +i 

136 A{ Ip1)=Tia*pEpAr ( n 
D0141J=1»N0C 
IPO=J*SS+i 

141 A( Ipo)=VALPP^ J)*YlDpA( J> 

JX=N0C+1 

JY=N0C+N05 

D0l42J=»jX,jY 

K0P=J*SS+1 

142 A{fCOP)=-(OCC+oCD*(SRi+SR2*THETAl> ) 
JX-JY+1 

JY=jY+N0S 

D0l43J=JX»jY 

KOP=J*SS+i 

143 A(KOP)=-(OCT+oCD»SR2) 

JX=JY+1 

JY=JY4-N0S 

D0144J=JX»jY 

KOPi=J*SS+i 

144 A{ KOP)=-(OCT+oCD) 

D081I=1»N0S 
D081J=1»N 
I0P=I*ZS+J*SS+1 
IF( ( J-I ) .NE.6)GoT 082 
A(IOP)=-1.0 
GOToSl 

82 ACIOP)=0.0 

81 continue 

Nl=NOS+l 

N2=2*NOS 

D083I=N1.N2 

D083J=1»N 

IOP=I*ZS+j*SS+l 

IF< ( I-J) ,NE.6)G0T084 

A( TOP) =-CONC 

GoTo83 

84 A(IOP)=0.0 

83 continue 



N1=N2+1 

N2=3*N0S 

D085I=N1,N2 

D085J=1»N 

IOP= I*ZS+j*SS+l 

IF( < I-J ) .NE.6 )G0To86 

A{Iop)=-CONT 

GOT085 

86 IF( ( J-I ) .NE.6)G0T087 
A( I0P)=-C0NT 
GoT085 

87 A(IOP)=0.0 
85 continue 

N1=N2+1 

N2=4*N0S 

DO 88 I=N1»N2 

DO 88 J=1»N 

KC= I*ZS+j*SS+l 

IF( j.GT.N0OGnT098 

lZ=l-36 

A(ICC)=COND*(Sp3-SR2 )*PI t IZ)/12.0 
GOTO 88 

98 IF ( (I-J) .NEoOi So fd 89 

aIkc)=-condi 

So To 88 

89 IFMI-J)*NE»18ISoTO90 
A(KC)=-COND2 
GOTO 88 

90 IF C<I-J).nE.6) Go To91 
a(ko=-cond 

GO To 88 

91 A{KC)=0.0 

88 continue 

N1=N24-1 

DO 92 J=l»NOC 

KC=Nl*ZS+j*SS+i 

92 A(!<C)=-1.0 
NS=N0C+1 

DO 93 J=NS,N 
KC=N1*ZS+J*SS+1 

93 A (KC )=0.0 
N1=N1+1 
Do 94 J = i 
KC=N1*ZS+J*SS+1 

94 A(lCC) = <AR2-AR3)*TApf 
NSi=NOC+i 
NK=NOC+NOS 

DO 95 J=NS,NK 
KC=N1»2S+j*SS+1 
P5 ATkC )=AR2*THETAi + ARi 



ny=nk+nos 

NK=NK+1 

DO 96 J=NK»NY 

KC=Nl-i^ZS+j*SS+l 

96 A(KC)=-(1.-AR?) 
nz=ny+nos 

NY=NY+1 

DO 97 J=NY,NZ 

KC=N1*ZS+J*SS+1 

97 A{KC)=-1,0 

N11=N1+1 
N12=N1+N0S 
DO 121 I=Nn tNia 
NlVsI-Nl 
DOl 26 J = 1 i N 
IAO=I*-ZS+j*SS+i 

if(j.gt,noc) go To 122 

A{ IAq) =-Wr ( j»NiY ) /X2»0 P I ( Ml Y ) -x-THETA^ / 1 ^ .<5 

GOTO 126 

122 IF( { I-.J) .nE. 44) 60 td 123 
AC iA0)=THETAi*THETA2 

Gotoi26 

123 IFC CI-J) ,mE.32) Go To 124 
AC IA0)=THETA2 

G0T0126 

124 A(lAO)=0,d 

126 CONTINUE 

121 continue 

N11=N12+1 

N12=N12+N0C 

D0127I=N11#N12 

D0128J=1»N 

IA0=I*ZS+J*SS+1 

IF( ( I-J).E0.62)G0T0129 

A { lAO) =0*0 

6oT0128 

129 A(IAO)=-1.0 

128 continue 

127 Continue 

CALL SIwplX{ZS»SS*m»mi*FALL»PP 1 »rPRl»PR2»IPP2» 
ILl* JL1»L2»JL2) 

PRlNTiOOO»A( 1) 

lOCO FoRMATC iX»*ToTAl benefit iis •*-»El ) 

D0181I=1*M1 
J=I + 1 

PRInT2000,pr2{ I) *A( J) 

181 continue 

2000 FoRmATC iX» l5 »* *»E16 ,a1 

stop 

END 
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IBFTC SImPLX 

IN THE program THE SIMPLEX TABLEAU IS STORED AS A SINGLE 
ARRAY MATRIX COLUMNWISE AND JHE OPTIMAL SOLUTION IS OB0A4-E+ 
IN THE SAME STORAGE LOCATION 

SUBROUTINE SIMPLX (ZS»SS»N, Ml , fall, PR3 , I PR1,PR2, I PR2» 
1L1»JL1>L2»JL2) 
commona ( 3OOO) 

dimension PRl ( IPRl ) »PR2t IPR2) ,L1 < JLl ) *1 2 ( JL2) 
integer ZS,SS,V,R,S,PR1 »PR 2 *FALL 
real' MAX 
R = 0 

Do 1 K= 1 »N 
Ll(K) = K 

1 PR1(K}=K 
LIO = N 

DO 2 I=1*M1 

2 L 2 ( n = I 
L20=M1 
D03I=liNi 

3 PR2il)sN+I 
GOT 01 O 3 

1 O 2 I S=PRl{Kpr 

PRl(KP)=PR2tlpi 

PR2nP)=S 

1 0 3 CALL MP7(0»ZS,SS,KP»L1,LiO,jli »MAX j 
IF(MAX.GT.0.3G0t014 

FALL=0 

return 

14 CALL MP 2 CL 2 M 20, JL 2 »IP»ZS,SS,KP»Q 1 »N»V) 

104 IF( Ip.NE.0)Gofol5 

FALL=: 1 '- 

return 

15 CALL MP3<0,M1»0,N,IP,Kp,zS,SS,1,U) 

GoTo IO 21 

END 


IBFTC MP 7 

subroutine MP7 determines the column that enters the basis 
subroutine MP7tZNR»ZS*SS,KP,Ll ,L10.JL1 ,MAX3 
ComMONA ( 3OOO) 

DIMEnSIONLI ( J tl i 
integer ZNR»ZS»SS 
real MAX 

kh=znr*zs+li < I )* ss+i 
MAX=A (KH ) 

KP=L1 «1 ) 

IF (L10.lt. 2) return 
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DOlK=2»LlO 
K H=ZnR*^S+l 1 ^K) *SS+i 
IF(A(kH).lE.MAX) goToi 
MAX=A(KH) 

KP=L1(K) 

1 continue 
return 
end 


IBFTC MP2 

subroutine MP2 DETERMINES THE ROW THAT LEAVES THE BASIS 
subroutine MP2(L2»L20,JL2»IP,ZS»SS»KP,Qi,n»IV) 

COMMONA (3OOO) 

DIMENSION L?IjL2j 
integer ZS,SS,Z 
V=IV 
IP=0 

IF«l20.lT.i ) return 

D0ll=l»L20 

tCH = L2 ( I )*ZS+i 

KH1=K,H+KP*5S 

1 IF (V*A( KHi ) .Gt.O* JGoT02 
RETURN 

2 Ql=V*A ( kH) /A< kHl ) 

IP=L2 <1 ) 

Z=I + 1 

if(z.gt.l20) return 
D 03t=Z»L20 
!CH=L2 ( n*zs+i 
khi=kh+kp»ss 

IF(V#A(KH1).LE.0,) 6oTo3 
Q=V*A(K:Hy/A(KHl ) 

IF(Q,GE.Q1) GoTOA 
IP=L2il) 

01 = Q 
GOT03 

A IF(Q.NE.oI) GoT03 
io=l 2 ( n 

D05K=1»N 

KHO=IP*ZS+K*Ss+i 

KH2=:IP»ZS+KP*SS+1 

KH=I0*ZS+k*SS+1 

qP=V*A(KH0)/A(KH2) 

qO=V*A(KH)/A(KH1) 

IF(QP.LT.qO) G0T03 
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5 IF(q0.lT.qp)GoT06 

6 IP = I0 

3 continue 

return 

END 


IBFTC MP3 

SUBROUTINE TRANSFORMS THE SIMPLEX TABLEAU 

subroutine MP3(I0»I1»K0,K.i »IP»KP,ZS,SS»P1»P2} 
COMMONAC 3000) 

Integer zs,ss,pi ,p2 

KH=IP*ZS+KP^SS+1 . 

PIV=1./A{«H) 

110=10+1 
111=11+1 
KK0=K0+1 
Ki<l = f<l + 1 

d6iii=ii0,ii1 

1=11-1 

IFd.EQ.IP) GOTOI 

KH0=I*ZS+KP*SS+1 

IF(P2.EQ.l) A(KHO) +A ( kHo ) »P I V 

Do2KK=Kk6»kIC1 

K = ICK-1 

IFifC.EO.KP) 6ot02 
KHl = I*ZS+K»SS4-i 
K.H2 =IP*ZS+K*Ss+1 
A(KH1)=A(KH1) - A(KH2)*A(kH0) 

2 continue 
1 continue 

IF(P1.nE,i) GoToA 
DOS KK=KK0,KKi 
K=KK-1 

kh2=ip*zs+k:*ss+i 

5 IF,(|C.nE*K:P ) A{tfH2) - -ATkM2)*pIV 

4 IF(P2.EQ,1) A(fCH)=PTV 

return 

END 
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APF3UDIZ - C 


A Discussion on the Results of the Cropping Policies 

Tile solutions reported for the cropping policies 
in Section 5-2 can he explained by the considerations of 
the annual return from each crop for every acre-feet of 
water supplied. Table 0-1 gives these informations. 

TABIS C - .1 


Crop. 

Annual 

water 

requirement 
in inche s 

Average””"' 
Annual 
yield is 
lb s/acre 

Value 

crops 

Es./lb 

of “ Annual 
in value 
in 

Ss./acre 

Annual 

lie turn per 

acre-feet 

0 f water 
supplied 

(6)=|-^x12*Q 

in lis./acre- 
feet 

(1) 

(2) 

(3) 

(4) 

(5) = 

(3' :-(4) 

“Wheat 

21 

900 

0.50 

450 

257 

Rice 

40 

800 

1.00 

800 

240 

Mai ze 

24 

1100 

1.00 

1100 

550 

Ootton 

29 

200 

0.50 

100 

42 

Sugar- 

50 

3000 

0.75 

2250 

540 

cane 






Other 

24 

600 

0.50 

300 

150 

crops 







¥e find from 

above table that 

'maize* and 



#■ t 


’sugar-cane' have the highest annual return per acre-feet 
of water supplied. They therefore enter the solution 



at their maximum' permissi'ble level, Npxt in order is 
wheat and we find that 'wheat' also enters the solution 
at its maximum permissible level. The 'other crops' 

i 

although it has a lov/er return than rice; also^ enters the 
solution at its maximum permissible level, because its 
seasonal water requirements are such that they can be 
met more or less by natural precipitation alone, which is 
available free. Water that is available after the 
requirements of the above crops are met, is made avail- 
able for rice and 'rice' enters the solution at a level 
lower than its maximum permissible level. 'Ootton' which 

I 

has a very poor return as compared to all the other crops 
considered does not enter the solution at all. 



